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Charmonium production in pp annihilation: 
Estimating cross sections from decay widths 
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The cross sections for the charmonium production processes pp — > rnW, where m is a light meson 
and $ is a charmonium state, are of great interest for the planned pp experiment PANDA. In 
this paper we estimate these cross sections using known results for the decays $ — » mpp, which are 
related to these reactions by crossing. In lieu of detailed experimental data on the decay Dalitz plots, 
we assume a constant amplitude as a first approximation; this implies a simple relation between 
the cross sections and decay widths. The single measured exclusive cross section of this type is 
pp — > -K°J/ip, which was reported by E760 to be 130 ± 25 pb near \fa = 3.5-3.6 GeV. In comparison, 
our constant amplitude estimate is about 300 pb at this energy. This suggests that our approach 
is useful as a simple estimate, but should be refined through detailed modeling of resonances and 
other energy-dependent effects in the experimental 9 — * mpp Dalitz plots. 

PACS numbers: 11.80.-m, 13.25.Gv, 13.75.Cs, 14.40.Gx 
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I. INTRODUCTION 

The proton-antiproton annihilation experiment 
PANDA, which is planned for the FAIR facility at GSI, 
is scheduled to begin data taking in the beginning of 
the next decade. PANDA is an internal fixed target 
experiment, with a design luminosity of 2 • 10 32 cm _2 s _1 
and a momentum range of Pp; t> = 1.5-15 GeV/c. This 
covers a mass range in which many interesting hadron 
resonances are anticipated, including glueballs, light 
hybrids, charmonia and charmonium hybrids pj. 

The physics goals of PANDA include a detailed investi- 
gation of the spectrum of charmonia and charmonium hy- 
brids, including determinations of masses, widths, quan- 
tum numbers and decay properties. Since the known 
charmonium spectrum is largely characterized by clear, 
well separated states, unlike the broad, overlapping res- 
onances of light quark hadron spectroscopy, it is antici- 
pated that this will be a relatively clean sector in which 
to search for higher-mass excitations such as hybrids and 
other unusual hadron resonances. 

Studies of conventional meson spectroscopy using 
proton-antiproton collisions have a considerable advan- 
tage over electron-positron colliders, in that all nonexotic 
quantum numbers can be accessed in formation reactions. 
Although J PC -exotics cannot be made in formation in pp 
annihilation, they can be made in associated production 
processes such as pp — > ttH c (where H c is a J PC -exotic 
charmonium hybrid). 



Production of the very characteristic J -exotic char- 
monium hybrids is one of the main goals of the PANDA 
project. These states are predicted to lie at rather high 
masses; the J PC = 1 h state, which is expected to be 
the lightest exotic charmonium hybrid, is predicted b 
the flux tube model HUSH and lattice Q CD Hf ' 
to have a mass in the range 4.2-4.4 GeV. This exotic 
resonance could be produced at PANDA together with 
a light meson m, such as an ri or ir°. (Both isospins are 
allowed in pp — ► m*!! because pp is a superposition of 1=0 
and 1=1.) 

Numerical estimates of charmonium and charmonium 
hybrid production cross sections are evidently crucial for 
the PANDA physics program, both to formulate detec- 
tion strategies and to evaluate luminosity requirements, 
as well as for detailed detector simulations with theoret- 
ically preferred final states. Unfortunately there are few 
estimates of these cross sections at present, and only one 
has been measured, pp — > ir°J/ip. (This result was re- 
ported by the E760 collaboration Q.) The purpose of 
this paper is to derive and evaluate simple estimates of 
these cross sections from measurements of charmonium 
partial widths to mpp three-body final states. 

The approach we follow in this paper is to use a mea- 
sured charmonium partial width into a three-body mpp 
final state, T^,^ mp p, to estimate the corresponding pro- 
duction cross section a p p^ m qi. The (presumably com- 
plicated) quark-gluon dynamics is subsumed in the un- 
known amplitude A4, which is common to both processes. 
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II. RELATING WIDTHS TO CROSS SECTIONS 



A. Three-body decay "if — > mpp 




FIG. 1: The three-body charmonium decay process ^ — > mpp. 



The three-body partial width of an initial charmonium 
state \1/ is an average over initial and sum over final po- 
larizations of a squared decay amplitude, multiplied by 
three-body phase space. The differential decay rate, in 
the form familiar as a Dalitz plot, is 
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1 



1 



mpp 



2S* + 1 (2tt) 3 32M| 



dm mp dm pp 



Here Sx is the spin of particle X, m 2 j = (pt + Pj) 2 is 
the invariant mass squared of the two-particle system 
ihj) [13 ■ an( i the squared amplitude in this formula is 
summed over all 4(2S* + l)(2S m + 1) initial and final 
polarization states. (This summation convention is con- 
venient for relating decays to cross sections.) 

As a simple approximation one can assume a constant 
decay amplitude A4, in which case the decay width is 
proportional to the area Ad of the Dalitz plot; 



1 



1 



>mpp 



where 



2SW + 1 (2tt) 3 32Af, 



A D = I I dm l mp dm 2 pp 



(2) 



(3) 



where t is the Mandelstam variable t — (p p — p m ) 2 ■ A 
factor of 1/4 has been included for the average over the 
four initial pp polarization states. 

The total cross section involves an integral over t be- 
tween the limits 



t±=m 2 p + m 2 m 



2(E p E m =F PpPm) C 



(5) 



Assuming a constant amplitude we can trivially evaluate 
this integral, with the result 



CT m) 



pp — ►m 1 !' 



l P 



647T p pt 



m cm — 
S 



(6) 



C. Connecting charmonium decay and production 



The partial width 



and the production cross 



section (r pp ^ m qi given above are simply related in the 
constant amplitude approximation, since both are given 
by simple kinematic and spin factors times the same spin- 
summed squared amplitude. Eliminating the common 
squared amplitude, we find the following relation between 
the cross section and decay width: 



2 M% 

&pp— >m* — 47T (2S , ip + 1) — rv[/_^ TO pp 

Ad 



Pm cm — 1 

S 

Pp cm 



(7) 



Note that in this approximation all the energy depen- 
dence in the cross section comes from the expression in 
square brackets. A similar approach has been used to 
estimate charmonium cross sections of relevance to QGP 
experiments |ll|. 



D. Experimental — > mpp decay widths 



B. Two-body production reaction pp — -> m\& 




FIG. 2: The related associated charmonium production pro- 
cess pp — > mty. 



In the unpolarized two-body production cross section 
the squared amplitude is again averaged over initial and 
summed over final states. The unpolarized differential 
cross section is given by 01 



da 
~dl 



256-7T \p pc 



^{ElMI 2 } (4) 



The estimate of the charmonium production cross sec- 
tions (7pp^ m * in Eq.(|7J) requires the three body decay 
partial widths T^,^ mpp . These branching fractions (and 
hence partial widths) have been measured for about 10 
mpp decay modes to date, all from J /ip or tp' initial states 
at e + e~ colliders. PDG values for the branching frac- 
tions are given in Table [I] Since the PDG compilation, 
new measurements of J /ip and ip' branching fractions to 
mpp have been reported by BESII [l^ and CLEO-c 0; 
these are summarized in Table ITTI 

Unfortunately, branching fractions of J PC ^ 1 — - char- 
monia to mpp final states have not yet been reported. 
Results for C = (+) states such as the r\ c and Xcj would 
be very interesting, as these may couple to mpp through 
a different mechanism than C = (— ) charmonia, per- 
haps with significantly larger amplitudes. This possibil- 
ity is suggested by the known four-body partial widths 
to TT + ir~pp, which are larger on average for the \J states 
than for the J/ip and ip'. 
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E. Previous Estimates 

The only previous estimate of near-threshold asso- 
ciated charmonium production cross sections we have 
found in the literature is the PCAC calculation of Gail- 
lard, Maiani and Petronzio 01- This reference con- 
siders the process ^ — > ir°pp, where $ is a generic 
charmonium state, and notes that the cross section is 
proportional to the peak on-resonance cross section for 
pp —> ^ — > all if one assumes that pion emission in the 
decay takes place from the p or p line in a two-step pro- 
cess, * — > pp — > TT°pp. Assumption of the usual PCAC 
coupling between nucleons and pions gives a relation be- 
tween the two cross sections, 



& pp— ►•n- * 



1 

4tt 



9a 
U 



j^tot max 

-L * &CP CTpp^^^allP^cm- 



(8) 



a CP is a dimensionless function of the cm. energy 
and the initial pp quantum numbers. Gaillard et al. 
estimated ucp numerically at a single pion energy, 
= 230 MeV, but did not evaluate it analytically. 
We have carried out this calculation using the angu- 
lar distributions given in their Eqs. (25-28). Defining 
Cm.n = (1/4) j_ 1 dcos(8) a m . n analogous to their Eq.(38) 
for acp, we find 
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(9) 



f(P), (io) 



/(/?) (11) 



where (3 = /3 w /3 p , 7 = 7^7^ and f(x) = taritT 1 (x)/x. 
The coupling ocp is a linear combination of these £ func- 
tions. For example, for a pure S-wave TT°J/tp system one 
requires a CP= — , 1 P± 1 H pp initial state, for which 
acp=- — £0,0 ■ For a CP= +, spin-triplet initial state, 
acp=+ is a linear combination of £1,0 and £1,1, 

a CP=+ = (1 - (S 3 2 » a, + (S 2 3 ) 6,1- (12) 

The spin matrix element {S 2 ) depends on the mix of 
L = J ± 1 waves in the initial pp state, and satisfies 
< (Si) < 1. (For an initial 3 Si pp state, (Sf ) = 2/3.) 

Numerical results for (Jpp^oj/^ using both this PCAC 
approach and the constant amplitude formula of Eq.Q 
are given below. 



III. NUMERICAL RESULTS FOR a pP ^ m9 
A. Method 

Using Eq.Q) and the partial widths given in Tables [I] 
and[Hl we have generated constant amplitude estimates 



for (jpp^m* for all cases with measured $ partial widths 
to mpp final states. The J/tp cross section estimates 
use PDG branching fractions and decay widths, and the 
ip' cross sections were calculated using BESII branching 
fractions for m = ir° and 77 and CLEO-c branching frac- 
tions for p, to and <f>. Our results, including combined 
errors from the total widths and branching fractions, are 
given in Table ITTTl 



B. The case pp — > w°J/ip 

The process pp — ► J/ i/jtt is of special interest because 
this is the single associated charmonium production re- 
action in pp annihilation for which we have an experi- 
mental measurement of the cross section. This reaction 
was studied by the E760 Collaboration at Fermilab [jj, 
who observed a few pp — > ir° <//?/> events as a background 
in their searches for 1 Pi h c and 2 3 Si rj' c charmonia. 
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FIG. 3: Kinematically allowed regions for the three-body de- 
cay J/ip — y TT°pp and the related charmonium production re- 
action pp — > ir°J/ip. 



The observed events correspond to a cross section of 
99 ± 40 pb just below 3.525 GeV and 156 ± 36 pb at 
3.61 GeV. This has also been quoted as a combined value 
of 130 ± 25 pb ^5|. The E760 measurements are shown 
in Fig0| together with the theoretical predictions. 

The constant amplitude cross section of Eq.J7|) is also 
shown in Fig0] This approximation predicts a = 299 pb 
at 3.52 GeV and 336p6 at 3.61 GeV, which overestimates 
the cross section at these energies by about a factor of 2- 
3. This suggests that the constant amplitude approxima- 
tion is useful as a first estimate, but should be improved 
through a more detailed description of the amplitude, for 
example by incorporating the contributions of individual 
resonances to the decay Dalitz plot. The importance of 
baryon resonances in this decay and in related \& — » mpp 
decays has been discussed previouslvin both experimen- 
tal [SH E3 and theoretical EE! references. 




FIG. 4: Theoretical and experimental cross sections for pp — > 
ir°J/tp. The theoretical predictions are the constant ampli- 
tude result Eq.lJJJ (solid) and the range of PC AC cross sec- 
tions, from Eq.@ (filled). The experimental points are from 
E760 0. 
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FIG. 5: Allowed regions for tp' - 
smaller physical region for J/tp 
ison (dashed, from FigHTl, 



7r°pp and pp — > TY°tp' . The 
Tv°pp is shown for compar- 



The PCAC result for the cross section for J/tp — > n pp, 
Eq.lJS}, is also shown in FigQ] Here we have assumed 
(rounded) PDG values of g A = 1.27 and fa = 0.131 GeV, 
PDG masses and widths, and a peak on-resonance cross 
section of cr r ^H J ^_ >aU = 5.0 ub |2Q|. Since the result de- 
pends on the initial pp quantum numbers assumed, we 
show the full set of values that follow from the range 
< (Sf) < 1 in Ea.iPJl. Evidently this PCAC pre- 
diction is also somewhat larger than the experimental 
E760 cross section in the 3.5-3.6 GeV region where we 
have data, and is rather similar to the constant ampli- 
tude result near threshold. The disagreement of PCAC 
with experiment may simply imply that other important 
processes interfere destructively with the nucleon pole 
diagram. 



It is notable that quite large maximum cross sections 
of 1-3 nb are predicted for ^ — J/tp and m = 77,0;, 7/. 
This suggests that it may be interesting for PANDA to 
consider reactions in which m is a light isoscalar meson. 
This is also supported by the relatively large branching 
fractions observed for charmonium decays to Tr + ir~pp. 
Although an enhancement of pp — > rrV& cross sections 
with light isoscalar m is an interesting possibility, one 
should note that the Dalitz plot for the decay J/tp — > rjpp 
(as one example) is crossed by N*(1535) bands (FigHJJ), 
and this baryon resonance contribution may be the rea- 
son for the large branching fraction. This is supported 
by results from BESII [HI. As these N*(1535) bands lie 
far outside the physical region for pp — > r/J/tp, it appears 
likely that the constant amplitude approximation Eq. J7J) 
will seriously overestimate this cross section. 



C. Other channels 

Our constant amplitude predictions for all pp — > 
cross sections are summarized in Table IIHI 

The pp — * mip' cross sections are typically predicted 
to be 10s of nb, at least an order of magnitude smaller 
than the corresponding J/tp cross sections. There is no 
indication of an enhancement of isoscalar m final states in 
this case. Although the tp' couplings are relatively weak, 
the decay tp' — > ir°pp (and the corresponding reaction 
pp — > 7r O- 0') may be the most interesting to study in 
the near future. As is evident in Fig|Sl the region of 
phase space sampled by this decay is much larger than 
in J/tp — y TT°pp, so more baryon resonance contributions 
can be investigated, and their effects on the extrapolation 
to the physical region for pp — > ir°tp' can be studied once 
data on this reaction becomes available. BESII 01 nas 
recently reported results on this tp' decay mode. 
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FIG. 6: Allowed regions for J/tp — > r\pp and pp — > rjj/tp, 
showing N*(1535) bands. (MzbF/2 limits are shown.) 
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IV. SUMMARY AND FUTURE 

We have used the partial decay widths of charmonia 
into mpp final states (where m is a light meson) and 
a constant amplitude approximation to estimate cross 
sections for the reactions pp — ► m^Sf, which are of great 
interest for the planned PANDA experiment. The sin- 
gle experimentally measured cross section of this type 
is pp — > TT°J/ip, which was reported by E760 to be 
130 ± 25 pb near 3.5-3.6 GeV. In comparison we predict 
a cross section of about 300 pb at this energy. Applica- 
tion of this constant amplitude approximation to other 
pp — * mJ/tp reactions suggests enhanced peak cross sec- 
tions of as large as 1-3 nb for light isoscalar m. This 
prediction may however be an artifact of large baryon 
resonance contributions to these specific J/tp — > mpp de- 
cays. These resonances need not contribute significantly 
to the production cross section if the decay resonance 
bands lie far from that region of phase space (see for ex- 
ample FigEI). The corresponding ip' decays do not show 
evidence of an enhancement of isosinglet m final states; 
the estimated pp — > mip' cross sections are 10s of nb. 

A next important step in these calculations will be to 
model the contributions of individual baryon resonances 
to the decays. In low-energy processes involving a sin- 
gle 7T°, such as J/tp — > ir°pp and pp — > ir°J/ip (Fig|3J), 
one may find dominance of the decay Dalitz plot by a 
few baryon resonances. In such cases, extrapolation of 
the amplitude to the reaction region may be relatively 
straightforward. Studies of the corresponding ip 1 pro- 
cesses may be instructive, since the region of phase space 
accessible to the decay is much larger (FiglSJl, and new 
data is expected from BESII and CLEO-c. In contrast, 
relating decays and production cross sections for heav- 
ier mesons such as the rj is clearly problematic, since the 
decay Dalitz plot provides a smaller window on the am- 



plitude, and the extrapolation from decay to production 
covers a much larger kinematic range (Fig[S|). For vector 
mesons such as the u> we encounter the additional com- 
plication that the couplings of baryon resonances to uN 
are not well established. The N* experimental program 
at TJNAF, which has a goal of determining these wN 
and related baryon resonance decay couplings, may also 
be useful in interpreting W — > mpp decays and clarifying 
their relation to charmonium production cross sections. 

Future measurements of mpp branching fractions of 
C = (+) charmonia such as the {xj} and r\ c will be es- 
pecially interesting, as these will allow estimates of their 
pp — > m^ production cross sections. As C = (+) char- 
monia employ different production mechanisms than the 
C = (— ) states, for example gg rather than ggg interme- 
diaries, they may have quite different production cross 
sections. 
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Decay 


~^mpp [^V] 


Experiment 


J/i> 




TT V pp 


99.2 ±8.9 


MRK2, MRK1, DASP 


J/i> 




VPP 


190 ± 18 


MRK2, MRK1, DASP 


J/i> 




p°pp 


< 28 


MRK2 


J/i> 




LOpp 


118 ±23 


MRK2, MRK1 


J/i> 




v'pp 


82 ±37 


MRK2, MRK1 


J/i> 




4>pp 


4.1 ± 1.4 


DM2 


4>' 




Tv'-'pp 


39 ±14 


MRK2 


V>' 




Ujpp 


22 ±9 


BES 


1>' 




<t>pp 


< 7.4 


BES 



TABLE I: Experimental charmonium partial widths to mpp fi- 
nal states. These numbers combine the PDG [Toj| total widths 
with the reported branching fractions. 



Decay 


r (BKB, [ e vi 

mpp L J 


plUlJSUjr y, 
mpp L v J 


,1 U - 
y) — > 7T pp 

1p — > T]pp 

if)' -> p°pp 

ip' — > upp 

■tp' — > 4>pp 


37.1 ±5.5 
16.3 ±3.8 


22 ± 12 
14 ±6 
17 ±8 
< 6.7 



TABLE II: — * mpp partial widths recently reported by 
BESII [p and CLEO-c [l|. These numbers combine the 
PDG 10] total width with the reported branching fractions. 



Reaction 




ESTIGeV] 


Ad [GeV 4 ] 


pp — > 7T U J/?/> 


420 ± 40 


4.28 


9.265 


pp — > J? J/V 


1520 ± 140 


4.57 


4.520 


pp-> p°J/il) 


< 450 


4.80 


2.114 


pp — > LJ J/?/) 

pp — > r/'J/ij) 


1900 ± 400 


4.80 


2.053 


3300 ± 1500 


4.99 


0.765 


pp — > J/V 1 


280 ± 90 


5.06 


0.452 


pp — > 7r u ?// 


55 ±8 


5.14 


30.500 


pp rjip' 


33 ±8 


5.38 


20.984 


Oil 

pp^ p i> 


38 ± 17 


5.59 


14.953 


pp — > LJ t/i' 


46 ±22 


5.60 


14.778 


PP — > 


< 28 


5.84 


9.118 



TABLE III: Charmonium cross sections from Eq.l|7|, given 
the ^ — > mpp partial widths of Tables [I] and [H] The size 
and energy of the cross section maximum and the area of the 
4' — > mpp Dalitz plot are given. These cross sections may 
be overestimated if the corresponding $ decay Dalitz plot is 
dominated by baryon resonances (see text). 



